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RESEARCH MEMORANDUM 


THE LONGITUDINAL STABILITY AND CONTROL CHARACTERISTICS OF 
A 60° DELTA-WING MISSILE HAVING HALF-DELTA TIP CONTROLS 
AS OBTAINED FROM A FREE-FLIGHT INVESTIGATION AT 
TRANSONIC AND SUPERSONIC SPEEDS 
By Martin T. Moul and Hal T. Baber, Jr. 

SUMMARY 

A rocket-propelled model of a cruciform 60° delta-wing missile 
configuration having half-delta tip controls has been flight-tested in 
the Mach number range 0.70 to I.65. Longitudinal stability derivatives, 
control effectiveness, and drag characteristics were obtained by use of 
the square-wave pulsed-control technique. 

The lift curves were smooth and linear within the angle-of-attack 
range (±10°) of the test. No unusual trends were noted in the stability 
derivatives with the exception of the damping- in-pitch derivative which, 
contrary to theory, had a minimum value at a Mach number of 1. 31. The 
lift-curve slope at supersonic speeds was less than the modified 
s lender-body theory results by from 7 ho 9 percent. The aerodynamic 
center shifted rearward 8-percent mean aerodynamic chord at transonic 
speeds. Pitching effectiveness of the tip controls was maintained 
throu gh the Mach number range but was reduced at supersonic speeds to 
about one-half the subsonic value. The steady-state angle of attack 
per unit control deflection at supersonic speeds was only 44 percent 
of the subsonic value. 


INTRODUCTION 


One type of control surface proposed for use on triangular wings 
at supersonic speeds is the half-delta tip control. Such a control on 
a 60° sweptback delta wing (ref. 1) has provided satisfactory rolling 
effectiveness and indicated good possibilities for aerodynamic balancing. 
Additional control- force and hinge-moment characteristics for a number 
of hinge-axis locations have been reported in reference 2. In order to 
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provide longitudinal-control data, an investigation of a cruciform 
missile having triangular wings and tip control surfaces was undertaken. 

From the flight test of a missile configuration having 60° delta 
wings and half-delta tip control surfaces, the longitudinal stability, 
control, and drag characteristics for a Mach number range of 0.70 to 
1.65 are presented. The results are compared with those of a 60° delta- 
wing airplane configuration having constant-chord full-span controls 
and with theory. 


SYMBOLS 


c wing chord, ft 

c" wing mean aerodynamic chord, ft 

S total wing area in one plane including body intercept, sq ft 

t wing thickness, inches; or time, seconds 

W weight, pounds 

Iy moment of inertia about Y-axis, slug- ft 2 

g acceleration due to gravity, ft/sec 2 

a angle of attack, deg 

5 control deflection, deg 

an/s normal accelerometer reading, g units 

a^/g longitudinal accelerometer reading, g units 

a-^/g transverse accelerometer reading, g units 

V velocity of model, ft/sec 

V c speed of sound in air, ft/sec 

M Mach number, V/V c 

p mass density of air, slugs/cu ft 

p coefficient of viscosity, slugs/ft-sec 
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R 

Reynolds number, pVc’/i-L 

q 

dynamic pressure, lb/sq ft; or pitching velocity, radians/sec 

p 

period, sec 

b 

exponential damping coefficient in e"^^, per sec 

C L 

/a n a Z \ W 

lift coefficient, [ — cos a + — sin a — 

\g e J qS 

C D 

( *1 a n . \ W 

drag coefficient, ( - — cos a + — sin a) — 

\ S g /qS 

Cm 

pitching-moment coe f f ic ient , 

Pitching moment about center of gravity 

qScT 


trim lift coefficient 

a trim 

trim angle of attack, deg 

^Anin 

minimum drag coefficient 

(VDW 

maximum lift-drag ratio 

dC D /clC L 2 

drag due to lift 

a/b 

amplitude ratio 

cp 

phase angle, deg 

CD 

angular forcing frequency, radians/sec 

a 

— i— — , radians /sec 
57.3 dt' 
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Derivatives : 


^Ct 

Ci^ = per deg; 


dc L 

% - §T’ deg 


Sc 

ClDa = So 21 ' per deg; 


dC 


Cm 6 = &T’ per deg 


dCm 

C m = =, per rad; 

q A 

° 2V 


bC, 


c itia = 7Z> P er rad 

^ ae 

3 2V 


MODEL AND APPARATUS 
Model Description 


A sketch of the model arrangement is presented in figure 1. The 
body consisted of a 7- inch-diameter cylindrical section and ogival nose 
and tail cones and had an over-all fineness ratio of 13.1. Triangular 
wings were mounted on the body in a cruciform arrangement. 

Wing and tip-control detail are presented in figure 2. The solid 
magnesium wings had a modified hexagonal airfoil section of constant 
thickness; the plan form was triangular, the leading edge being swept 
back 59° 32*. The tips of two opposite wing panels were replaced by 
steel, all-movable, half-delta, control surfaces mounted on steel 
torque rods. The controls had a double— wedge section and a constant 
thickness ratio of 3 percent. There was a parting-line gap between 
the control and wing of 0.03 inch. The ratio of exposed wing area 
(inclusive of controls) to control area was 9:1- Photographs of the 
model and the wing-control arrangement are presented in figure 3- 

A hydraulic pulsing system supplied from an accumulator was uti- 
lized to provide a programmed square-wave deflection of the controls. 

The two control deflections were ±9.6°. At a Mach number of about 1.0, 
the pulse frequency was decreased by means of a switch, sensitive to 
total-head pressure, which decreased the voltage field of the programming 
motor. 

The model was boosted to supersonic velocities by means of a solid- 
propellant rocket motor of 6000 -pound thrust and 3-second duration. 

The launching technique was the same as that reported in reference 3. 
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The physical characteristics of the model are presented in the 
following tat>le: 

Weight and balance: 

Weight, lb • • • • 

Center of gravity, percent c back of leading edge of c . . . . 12.6 

I y , slug- ft 2 8 - 1 

Wing: o Ak 

Total wing area in one plane, sq ft 2.04 


Mean aerodynamic chord, ft 

Thickness ratio at wing-body juncture 0.03 

Tip control: 

Area in one plane, sq ft 0.19 

Mean aerodynamic chord, ft 0.3o 


INSTRUMENTATION 


The model was equipped with an NACA 8-channel telemeter which 
transmitted continuous records of normal (2 ranges), longitudinal, and 
transverse accelerations, angle of attack, control deflection, total 
pressure, and a calibrated static pressure. Angle of attack was measured 
by a free-floating vane extended from the nose on a sting and total 
pressure by a total -pressure tube extended below the body. The posi- 
tions of these two instruments are shown in figure 1. 

The model trajectory was determined by a modified SCR-58^ type 
radar tracking unit. Model velocity was obtained from a CW Doppler 
velocimeter and from total and static pressures. A radiosonde released 
at the time of flight measured temperatures and atmospheric pressures 
through the altitude range traversed by the model. 


TECHNIQUE 


As the model decelerated through the Mach number range, longitudi- 
nal responses were produced by a programmed square-wave deflection of 
the tip controls. The angle of attack measured at the nose was cor- 
rected to the center of gravity by the methods of reference 4. The 
normal acceleration and corrected angle-of-attack oscillations were 
analyzed by the methods of reference 3 to determine the stability and 
control derivatives. 
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When the model drag was analyzed, the lift range at supersonic 
speeds was not large enough to yield experimental values of maximum 
lift-drag ratio. At these Mach numbers calculated values of (L/D) max 

were obtained from 


1 



which assumes the extensions of the lift-drag polars to vary as a 
parabola. 

Although the model encountered transverse accelerations up to 
±2.5g during the flight, no effects of yaw were noted on the longitudi- 
nal responses. 


ACCURACY 


The measured quantities are believed to be accurate within the 
following limits: 


M 


Limit of accuracy of 


M 

a 

6 

c L 

C Drain 

0.8 

±0.07 

±0.4 

±0.2 

±0.019 

±0.010 

1.6 

±.03 

±.4 

±.2 

±.oo4 

±.002 


These errors, dependent upon telemeter and radar precision, are 
primarily systematic in nature; parameters which are dependent upon 
differences in measured quantities, such as C^ and Aa/A 5 , are more 

accurately determined than the above errors indicate. 
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RESULTS AND DISCUSSION 


The scale of the flight test is presented in figure h; Reynolds 
number, based on wing mean aerodynamic chord, varied from 5*5 X 10^ 

to 16.3 X 10^. Sample time histories of the model flight at subsonic 
and supersonic speeds are presented in figure 5» The Mach number vari- 
ation with time, programmed control deflection, measured angle of 
attack, normal, and longitudinal accelerations are presented. Differ- 
ences in pulse duration, model natural frequency, and angle-of-attack 
range between subsonic and supersonic speeds may be noted. 


Static Stability 

Pitching moment .- In figure 6 is presented the variation with 
Mach number of the period of oscillation. From 0. 602 second at 
M = 0.70 the period decreased to 0.205 second at M = I. 63 . 

The static stability derivative Cn^ was reduced from the period 
of oscillation and presented in figure 7- A transonic increase in C^, 

of about -0.005, is noted between Mach numbers 0 . 90 and 1.00. At super- 
sonic speeds decreased gradually from -O.OI 36 at M •= 1.0 to 

-0.0083 at M = 1.63. 

Lift . - From lift-coefficient and angle-of-attack time histories, 
lift-curve plots were constructed for determining lift-curve slope Cp^. 

Some typical plots are presented in figure 8 for control deflections 
of ±9*6°. The lift curves are smooth and linear within the angle-of- 
attack range of this test. 

The slopes of the lift curves were measured and are presented in 
figure 9. Lift-curve slope increased from 0.0^38 at M = 0.72 to a 
peak value of 0.0525 at M = 1.02. As the Mach number continued to 
increase, decreased to O.O 387 at M = 1 . 62 . 

At a Mach number of 1. 62 a comparison is made with unpublished 
data from tests conducted in the Langley 9- inch supersonic tunnel of a 
similar wing-body configuration differing in body length. The wind- 
tunnel configuration was longer by a scaled 25 - inch cylindrical section 
inserted between the wing and the nose. Good agreement is shown in 
the Ci^ of these models. 

The lift-curve slope of a 60° delta-wing airplane reported in 
reference 5 is presented for comparison. Comparable body-upwash effects 
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are expected for the airplane and the missile since the body-diameter — 
wing-span ratios were nearly equal, being 0.226 for the airplane and 0.219 
for the missile. The lift-curve slopes agreed favorably at supersonic 
speeds, but a difference is noted at transonic speeds. This difference is 
attributable to the critical Mach numbers, the airplane having a lower 
critical Mach number because of its 65-series airfoil, 6.5-percent-thick 
section. 

In reference 6 Nielsen and -Kaattari have modified the s lender- 
body-theory results for lift-curve slopes of wing-body combinations 
having cylindrical afterbodies. From figure 9, theory is seen to pre- 
dict the lift-curve-slope trend, with the test values being 7 to 9 per- 
cent lower than the theoretical values at all Mach numbers. This dif- 
ference is believed attributable to the basic-wing experimental lift- 
curve slope being less than theory (for example, ref. 7) rather than to 
any negative lift realized by the tail cone. 

Aerodynamic center . - The aerodynamic-center location in percent 
of c is presented in figure 10. Aerodynamic center was at 31 percent c 
at subsonic speeds and shifted rearward in the transonic range to 
39 percent c at M = 1.05- With further increasing Mach number, the 
aerodynamic center shifted gradually forward. The transonic aerodynamic- 
center shift was only 8 percent c for this configuration as compared 
to 12 percent c" for the canard configuration of reference 3- In 
reference 7 the aerodynamic center of a wing-body combination having a 
cylindrical afterbody is determined from modified s lender-body theory. 

Good agreement is shown between this test and theory and indicates that 
the loss of lift on the tail cone is slight as suggested in the lift 
discussion. Also presented are the aerodynamic-center locations for a 
delta-wing airplane having a shorter forebody (ref. 5) and a delta-wing 
missile having a longer forebody. As expected, the configuration with 
the shortest forebody had the most rearward aerodynamic-center location. 
The change in missile forebody length caused an aerodynamic-center 
movement of 15 percent c" as compared with a predicted movement of 
10 percent c* from the modified s lender-body theory. 


Dynamic Stability 


The exponential damping constant b is determined from the ampli- 
tude ratios of the decaying oscillations by the methods of reference 3 
and is presented as a function of Mach number in figure 11. The expo- 
nential damping constant varied from 1.10 at M = 0.70 to 3-0 at 
M = 1.55 with a slight decrease occurring between M = 1.10 and 1.30. 


The 


The damping- in-pitch derivative C m ^ + is shown in figure 12. 

trend of C™ + C^. with Mach number with minimum points occurring 
m q nL a 
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at Mach numbers 0.97 and 1.33, is unusual. With increas ing Mach number 
above 1.33 C m + increases rapidly to a maximum value of -3-33 

at a Mach number of 1.62. The trend through the Mach number range 0.75 
to 1.10 is similar to that of the canard missile of reference 3 in that 
the minimum and maximum points occur in the same order. Modified 
s lender-body- theory results of reference 8 for a delta-wing— body con- 
figuration with no afterbody predict approximate magnitudes but do not 
predict the Mach number trend. The damping- in-pitch derivative of a 
60° delta-wing airplane configuration (ref. 5) is seen to be greater 
than the missile damping through most of the supersonic Mach number 
range. From a consideration of center-of- gravity locations (20.7 per- 
cent ~ for airplane and 12.6 percent c for missile) and body lengths, 
a greater damping- in-pitch derivative would have been predicted for the 
missile than for the airplane. The accuracy of % * is question- 

able since the lift and damping- in-pitch terms which contribute to the 
total damping (ref. 3) are approximately equal in magnitude. Any error 
occurring in C]^ or b will cause a correspondingly larger error in 

Cm q + 

It should he noted that although the damping- in-pitch derivative 
of wing-body combinations may be extremely small compared with that of 
some present-day wing- tail missile configurations, the differences in 
total damping are not so great. For example, a comparison of the tail- 
less model of this paper with the canard missile configuration of refer- 
ence 3 shows that, although the ratio of damping- in-pitch derivatives 
is 10:1, the ratio of total damping is about 3:1* 


Control Effectiveness 

Longitudinal trim curves .- The trim angles of attack are presented 
in figure 13 for the two programmed control deflections. For the con- 
trol deflection of -9.6°, a^rim decreases from 4° at M = 0.7 to 
2.1° at M = 1.15 and then increases to 2.6° at M = 1.60. The larger 
trim angles of attack at subsonic speeds result from the smaller static 
stability at these speeds. The differences in a trim ^ or the two e< l ua l 
control deflections indicate an out-of-trim condition which is possibly 
due to slight asymmetries in the model and in the angle-of-attack 
indicator. 

In figure lb is presented the trim-lift coefficients for the two 
control deflections. The C]^ . curves resemble the o^trim curves, 

the values decreasing with increasing Mach number. There is an out-of- 
trim lift coefficient as there was an out-of-trim angle of attack. 
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The trim lift per trim angle of attack of this configuration is 
about 30 percent less than that of the canard configuration because of 
the negative C ra& of the tip controls. 

Trim angle of attack and lift per unit control deflection . - From 
the ^t r im curves the trim angle of attack per unit control deflection 

is determined and presented as a function of Mach number in figure 15. 
The Ax-h-r im/ AS decreases abruptly near a Mach number of 1.0 and, at 

supersonic speeds, is only 44 percent of the subsonic value. 

The change in trim-lift coefficient due to a unit control deflec- 
tion is presented in figure 1 6. From 0. 0166 at M = 0. 7, the value of 
ACL t r in] /A5 decreases rapidly to O.OO65 at M = 1.40 primarily because 

of the rearward aerodynamic-center shift at transonic speeds. 

Control lift and pitching effectiveness .- The control lift effec- 
tiveness Cx,g including the lift induced on the wing was determined 

from the trim conditions of angle of attack and lift coefficient as in 
reference 3 and is presented in figure 17. From 0.0039 at M = 0.70, 

Cl 5 increases with Mach number to a peak of 0.0062 at M = 0. 91 and 

then decreases gradually to 0.0022 at M = 1.62. Lagerstrom (ref. 9) 
has obtained from the linearized theory the lift of a half-delta tip- 
control surface on a triangular wing and the corresponding wing- induced 
lift. The linearized theory does not predict the experimental results 
although the Mach number trends are similar. 

Control-surface pitching effectiveness is presented in fig- 

ure 18. Pitching effectiveness remains nearly constant at subsonic 
speeds but decreases gradually at supersonic speeds. At a Mach number 
of 1.5, C m g is about 56 percent of the subsonic value. 

The Cm^ of a 60° delta-wing airplane having constant-chord, full- 

span elevons (ref. 5) is also presented, after being corrected to the 
center-of-gravity location for comparable static stability. The ratio 
of control area to total wing area for the airplane is 2.46 that of the 
missile so that a larger C m g, based on total wing area, is expected 

for the airplane. The ratio of C m5 of the airplane to C m& of the 

missile varies from 3.6 at a Mach number of 0.9 to 2. 6 at a Mach number 
of 1.6 and indicates that the trailing-edge flap is the more effective 
control. Over the Mach number range, however, the loss of effectiveness 
is greater for the trailing-edge flap and at M = 1.6 the pitching 
effectiveness of the two controls is about equal. 
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Drag 


The variation of total drag of the model with lift coefficient is 
shown in the drag polars of figure 19 for a Mach number range of 0. 75 
to 1.55. The subsonic and supersonic polars indicate that the induced 
drag varied parabolically with the lift; however, this variation was 
altered slightly in the transonic range due to the rapid increase in 
wave drag with Mach number. 

Minimum drag coefficient as a function of Mach number is shown in 
figure 20. This is the drag at zero lift of the model with controls 
at 9*6°. A sharp drag rise occurred at transonic speeds reaching a 
peak of = 0.055 a t M = 1.10. The minimum drag coefficient 

decreased gradually through the supersonic region to a value of 0.0^3 
at M = I.55. 

The variation of drag due to lift with Mach number is shown in 
figure 21. A subsonic value of 0.55 was obtained for dCp/dCL 2 at 

M = 0.7. The trend showed a decrease to 0.^3 at M = 1.0 and a slight 
increase for the tested supersonic range. The same trend was indicated 
for 1/57*301^ with a slightly greater increase at supersonic speeds. 

The maximum lift-drag ratio, shown in figure 22, was measured 
directly from drag polars at subsonic and transonic speeds. At super- 
sonic speeds, the condition of maximum lift-drag ratio was not attained 
by the model. Calculated values of (L/D) max were obtained for these 


variation of drag with lift. The data indicate that (L/D) max decreased 

from 6. 3 at M = 0.70 to 3-2 at M = 1.1 and then gradually increased 
to 3.6 at M = I.55. 


Comparison of the results of two methods for determining longitudi- 
nal frequency-response characteristics from transient-response data is 
made in figure 23 for two Mach numbers. One method employs the experi- 
mental stability derivatives and the linearized equations of motion. 

The second method, as outlined in reference 10, utilizes the Fourier 
synthesizer to analyze the experimental transient-response data. In 
the analysis, amplitudes at 2k equal time increments were used to 
represent the transient response. 

For both Mach numbers the resonant frequencies as determined by 
the two methods are in good agreement and indicate the linearity of the 
static stability derivative. The peak amplitudes are a function of the 


Mach numbers by using 



assuming a parabolic 


Frequency Response 
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system damping and in figure 23 indicate a difference in damping. From 
the amplitudes the damping ratios (percent of critical damping) were 
obtained and were found to be only slightly different for the two 
methods. For M = 1. 60 the damping ratios were 0.11 and 0 . 13 , respec- 
tively, for the calculated and synthesizer methods. The same ratios 
were 0.10 and 0.13 tor the M = 0 . 7 ^- case. At M = 1. 60 the phase- 
angle curves are noted to diverge at higher frequencies. For a linear 
system of this order, the phase angle should approach -180° at these 
frequencies (verified by the calculated results). The trend of the 
synthesizer phase-angle curve suggests the presence of higher order 
terms in the transient response or a limitation of the method to a 
lower frequency range. 


Eesults from a flight investigation at Mach numbers 0. 70 to 1.65 
of a 60 ° delta-wing missile configuration having half-delta tip controls 
indicate the following conclusions: 

1. No unusual trends were noted in the stability derivatives with 
the exception of the damping- in-pitch derivative which, contrary to 

s lender-body theory, decreased to a minimum value at a Mach number of 
1.33 and then increased with increasing Mach number. The lift curves 
were smooth and linear within the angle-of-attack range (±10°) of the 
test. The lift-curve slope was in good agreement with that of a tail- 
less airplane configuration but was less than the modified s lender-body- 
theory results by from 7 to 9 percent. The aerodynamic center was 
located at 31 percent c" at subsonic speeds and shifted rearward 8 per- 
cent c at transonic speeds. At supersonic speeds the aerodynamic- 
center location was in good agreement with modified s lender-body- theory 
results . 

2. Pitching effectiveness of the tip controls was maintained 
through the Mach number range but the value at a Mach number of 1.50 
is only 5 6 percent of the subsonic value. Per unit area the trailing- 
edge flap was a more effective longitudinal control than the half-delta 
tip control. The steady-state change in angle of attack per unit con- 
trol deflection Aoc/A5 at supersonic speeds was only bb percent of the 
subsonic value. 

3. The minimum drag coefficient was a maximum of 0.055 at a Mach 


CONCLUSIONS 


number of 1.10. The variation of drag with lift 



was larger 
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than l/57. 30]^ (where is the lift-curve slope) at all Mach num- 

bers. Maximum lift-drag ratio was 6.3 at subsonic speeds and about 
3.1 at supersonic speeds. 

c <■ >j£ J* 1 00X e V o 0,1 "iu 2 #Dl.no 8 njBTT . o rroscfuS 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 

Langley Field, Va. 
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Figure 1.- Model arrangement. All dimensions are in inches. 
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Figure 2.- Wing and tip-control detail. All dimensions are in inches. 
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(b) Wing and control- surface arrangement. 
Figure 3 »- Model photographs. 
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Figure 4.- Scale of flight test, based, on wing mean aerodynamic chord 

of 1.46 feet. 
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Time , /, sec 
(a) Supersonic. 

Figure 5.- Sample time histories of model flig 1 ' -1 ' 
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(b) Subsonic. 
Figure 5.- Concluded. 
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Figure 6.- Variation of period with Mach number. 
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Figure 7.- Variation of static pitching-moment derivative with 

Mach number. 
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(a) 5 = -9.6°. 



(b) 6 = 9.6°. 

Figure 8.- Typical lift-coefficient variation with angle of attack. 
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Figure 9.- Variation of lift-curve slope vith Mach number. 
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Figure 10.- Variation of aerodynamic-center position vith Mach number. 
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Figure 11.- Variation of exponential damping constant b with Mach number. 
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Figure 12.- Variation of aerodynamic damping- in-pitch derivative wit’ 

Mach number. 
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Figure 13.- Variation of trim angle of attack vith Mach number. 



Figure 14.- Variation of trim lift coefficient with Mach number. 
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Figure 15 . - Trim angle of attack produced by a unit control deflection. 
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Figure 16.- Variation with Mach number of trim lift coefficient due to 

control deflection. 
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Figure 17.- Control- surface lifting effectiveness. 



Figure l8.- Control-surface pitching effectiveness. 
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Figure 19.- Lift-drag pclars. 
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Figure 19.- Continued. 


NACA RM L52H14 



-4 ~3 :2 -/ O ./ .2 .3 .4 



-4 -3 -.2 -./ O ./ .3 .4 


Cl 


Figure 19.- Concluded. 
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Figure 20.- Variation of minimum drag coefficient with Mach number. 
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Figure 21.- Effect of lift on drag. 
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Figure 22.- Variation of maximum lift-drag ratio -with Mach number. 
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Figure 23.- Comparison of frequency-response methods. 
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Figure 23.- Concluded. 
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